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O . ABSTRACT 
Oh. 

\ We present a detailed study of molecular gas in the central 10 pc of the 

Galaxy through spectral line observations of four rotation inversion transitions 
of NH3 made with the Very Large Array. Updated line widths and NH3(1,1) 
opacities are presented, and temperatures, column densities, and masses are de- 
rived for the major molecular features. We examine the impact of Sgr A East 
^ on molecular material at the Galactic center and find that there is no evidence 

<N ■ that the expansion of this shell has moved a significant amount of the 50 km s 

' giant molecular cloud. The western streamer, however, shows strong indications 

that it is composed of material swept-up by the expansion of Sgr A East. Us- 
ing the mass and kinematics of the western streamer, we calculate an energy of 



O 



Qh' -E^at = (2 — 9) X 10^^ ergs for the progenitor explosion and conclude that Sgr A 

O ' East was most likely produced by a single supernova. The temperature struc- 

' ture of molecular gas in the central ~ 20 pc is also analyzed in detail. We find 

■ that molecular gas has a "two-temperature" structure similar to that measured 
by Hiittemeister et al. (1993a) on larger scales. The largest observed line ra- 

^ ' tios, however, cannot be understood in terms of a two-temperature model, and 

■ most likely result from absorption of NH3(3,3) emission by cool surface layers 
of clouds. By comparing the observed NH3 (6,6)-to-(3,3) line ratios, we disen- 
tangle three distinct molecular features within a projected distance of 2 pc from 
Sgr A*. Gas associated with the highest line ratios shows kinematic signatures 
of both rotation and expansion. The southern streamer shows no significant ve- 
locity gradients and does not appear to be directly associated with either the 
circumnuclear disk or the nucleus. The paper concludes with a discussion of the 
line-of-sight arrangement of the main features in the central 10 pc. 
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1. Introduction 

Since its discovery by Balick & Brown (1974), the bright (~ 1 Jy), compact radio source 
called Sagittarius A* (Sgr A*) has been hypothesized to mark the nucleus of the Milky Way. 
Today, the dynamical center of the Galaxy has been precisely determined through infrared 
observations of stellar orbits in the central 1" (0.04 pc) (Schodel et al. 2002; Ghez et al. 2003). 
The inferred mass of ~ 4 x 10® Mq contained within a radius of ~ 60 AU is some of the 
most convincing evidence for the existence of a supermassive black hole. Based on Very Long 
Baseline Array observations of the proper motion of Sgr A*, Reid & Brunthaler (2004) find 
that the peculiar motion of Sgr A* relative to the dynamical center of the Galaxy is — 0.4±0.9 
km s~^. Such a small velocity so near a supermassive black hole implies that Sgr A* must 
contain at least 10® Mq (Reid & Brunthaler 2004). It is therefore most hkely that that the 
radio emission from Sgr A* originates in either the inner region of a radiatively inefficient 
accretion flow (Quataert 2003) or the nozzle of a jet emanating from the supermassive black 
hole (Yuan et al. 2002). 

With a bolometric luminosity of ~ 10^® ergs s""^ (Narayan et al. 1998, and references 
therein), the emission from Sgr A* is extremely sub-Eddington {LEdd — 5 x 10^^ ergs 
for a 4 X 10® Mq black hole). However, the close proximity of Sgr A* (d — 8.0 ± 0.5 kpc, 
Reid (1993)), makes it uniquely suited for the detailed study of the environment around 
a supermassive black hole. In the central ~ 2 pc of the Galaxy, Sgr A* is surrounded 
by multiple arcs of ionized gas, which together are called the "mini-spiral" or Sgr A West 
(Ekers et al. 1983; Lo & Claussen 1983; Roberts & Goss 1993). These arcs are surrounded 
by an apparent "ring" of dense (~ 10^ cm~^) molecular material called the circumnuclear 
disk (CND, Glisten et al. (1987); Wright et al. (2001)). Neither of these features appears 
to be gravitationally stable, and it is estimated that (barring other external pressures) this 
material will accrete into the nucleus on the order of 10^ — 10^ yr (Wright et al. 2001). 

For the past two decades, the origin of the dense gas and dust in the CND and Sgr 
A West has remained unclear. Many attempts have been made to detect kinematic con- 
nections between two nearby giant molecular clouds (GMCs) and the CND. This work has 
focused primarily on observations of rotation inversion transitions of NH3 at frequencies of 
23-25 GHz. These transitions trace densities similar to those traced by the 3 mm HCN(l-O) 
transition inn,^ ~ lO^cm^'^), but the higher energies and lower opacities of the NH3 transi- 
tions significantly reduce the effects of self- absorption and absorption along the line-of-sight 
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(LOS) (McGary et al. 2001). Due to the close frequency spacing of these spectral hues, ob- 
servations of multiple transitions can be made with the same instrument and receiver making 
calculation of line ratios (and therefore rotational temperatures) relatively straightforward. 
In addition, hyperfine splitting of the NH3 lines enables direct calculation of gas opacity (Ho 
& Townes 1983). 

Previous NH3 observations detect a long filamentary "southern streamer" extending 
northwards from the "20 km s'^" GMC (M-0. 13-0.08; Giistcn ct al. (1981)) towards the 
the southeastern edge of the CND (Okumura et al. 1989; Ho et al. 1991; Coil & Ho 1999, 
2000). (Sec Figure Ic for positions of the main molecular features in the region.) Increased 
line widths and evidence of heating as the southern streamer approaches the Galactic cen- 
ter indicate that gas may be flowing from the 20 km s"^ GMC towards the circumnuclear 
region. Supporting morphological evidence for this connection has come from observations 
of HCN(3-2) (Marshall et al. 1995), ^^C0(2-l) (Zylka et al. 1990), and 1 mm continuum 
emission (Dent et al. 1993; Zylka et al. 1998). However, the southern streamer shows neither 
the velocity gradient nor curved morphology that would be expected for a cloud infalling 
towards the nucleus (Coil & Ho 1999). NH3 emission also becomes very faint near the CND 
making it difficult to study the gas in this region (McGary et al. 2001). Two additional 
connections originating in the northern fork of the 20 km s~^ GMC (Coil & Ho 1999, 2000) 
and the northern ridge (McGary et al. 2001) have also been suggested based on NH3 obser- 
vations. A fourth possible connection originating in the "50 km s~^" GMC (M-0. 03-0. 07) 
was detected in HCN(l-O), but it has not been observed in NH3 (Ho 1995). 

In this paper, we combine data from previous interferometric observations of four rota- 
tion inversion transitions of NH3 to study the molecular environment in this region in detail. 
The observations and velocity integrated maps, which were previously published in McGary 
et al. (2001) and Herrnstein & Ho (2002), are briefly reviewed in §2. We then derive physical 
parameters for the molecular clouds in §3, and use these results to investigate the impact of 
Sgr A East on molecular gas in the central 10 pc in §4. Based on the interaction of Sgr A 
East with the western streamer, we place limits on the energy of the progenitor explosion. 
In §5, we discuss the detailed temperature structure of molecular clouds at the Galactic 
center. Measured temperatures for clouds more than 2 pc from Sgr A* are consistent with 
a two-temperature structure in which roughly one quarter of the gas is contained in a hot 
(~ 200 K) component. This structure is similar to the model used by Hiittemeister et al. 
(1993a) to characterize temperatures on larger scales. NH3 (6,6)-to-(3,3) line ratios in the 
central 2 pc, however, exceed theoretical hmits for gas in LTE (Herrnstein & Ho 2002). In 
§5.2, we model these line ratios as resulting from absorption of NH3(3,3) by cool material in 
a shielded layer of the cloud. In §6, we disentangle emission from three distinct molecular 
clouds within a projected distance of 2 pc from Sgr A* through a comparison of NH3(3,3) and 
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(6,6) emission. We conclude this paper with a model of the hne-of-sight (LOS) distribution 
of the main features in the central 10 pc based on our molecular data and in the context of 
existing data at other wavelengths. 

2. Observations and Velocity Integrated Maps 

The large range of velocities observed at the Galactic center and the relatively narrow 
velocity coverage of current interferometers at centimeter wavelengths make spectral line 
observations of the region difficult. A velocity coverage of ±110 km s~^ is necessary to 
detect all of the emission from the CND, but additional clouds with velocities as high as 
-185 km s~^ have been detected near Sgr A* (Zhao et al. 1995). Using the Very Large Array 
(VLA)\ we observed the central 10 pc of the Galaxy in NH3(1,1), (2,2), and (3,3) in 1999 
March. The five pointing mosaic fully samples the central 5' (10 pc) of the Galaxy with a 
velocity coverage of of -140 to +130 km s~^. The resulting data are more complete in both 
velocity coverage and spatial sampling than previous interferometric observations of NH3 
emission in this region, which tended to concentrate on velocities associated with the GMCs 
(Coil & Ho 1999). However, the wide velocity coverage of our data is gained at the expense 
of channel width, which is only 9.8 km s~^. 

Figure la-c shows velocity integrated maps of NH3(1,1), (2,2), and (3,3) emission from 
the central 10 pc (see McGary et al. (2001) for details on the observational setup and data 
reduction). The position of Sgr A* (aaooo = 17M5"^40^0, ^2000 = -29°00'26".6) is marked by 
a star in each panel. Major molecular features discussed in this paper are labeled in Figure 
Ic. After application of a Gaussian taper to de-emphasize long baselines, the synthesized 
beam in all three maps is roughly 16" x 14" with a position angle of ~ 0°. The rms noise 
level, ajK, for each {J,K) rotation inversion transition is an = 0.28 Jy beam~^ km s~^, 
(722 = 0.30 Jy beam~^ km s"-*^, and (J33 = 0.33 Jy beam~^ km s~^ (McGary et al. 2001). 
Contours in all four panels of Figure 1 are in steps of 1.32 Jy Beam~^ km s~^, thereby 
facilitating comparisons between different transitions. However, all emission exceeding 3a jk 
is considered significant (see McGary et al. (2001) for maps showing 3a jk contours). In 
order to have a constant noise level across the map, we have not corrected for primary beam 
attenuation at the edge of our mosaic, and the data are not sensitive to emission ^ 150" 
from Sgr A* (McGary et al. 2001). The 50 km s^^ GMC, molecular ridge (composed of SEl 
and SE2), and southern streamer all extend past the edge of the mosaic (see Armstrong & 
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Barrett 1985; Dent et al. 1993; Zylka et al. 1998, for large-scale maps showing the extent of 
these clouds). However, the outer edges of both the northern ridge and western streamer 
occur well within our mosaic and reflect intrinsic cloud edges (McGary et al. 2001). 

As reported in McGary et al. (2001), we find dense molecular gas throughout much of 
the region. A notable exception is the far western and northwestern parts of the mosaic. 

These regions also lack emission from HCN(l-O) (Wright et al. 2001) and thermal dust (Dent 
et al. 1993; Zylka et al. 1998; Pierce-Price et al. 2000). Emission from NH3 is predominately 
associated with the 20 and 50 km s^^ GMCs. Both the southern streamer and northern fork 
(Coil & Ho 1999) arc extensions of the 20 km GMC. The two southeastern clouds (SEl 
and SE2) form part of the "molecular ridge" that connects this GMC to the 50 km s^^ GMC 
in the northeast (Coil & Ho 2000). Only the CND, northern ridge, and western streamer 
show velocities outside the ~ 20 — 50 km s~^ range, and are not direct extensions of these 
GMCs. Overall, we find some evidence for kinematic connections between giant molecular 
clouds in the central 10 pc, but the expansion of the Sgr A East shell appears to be a 
dominating force in the region, possibly pushing much of the molecular gas away from the 
nucleus (McGary et al. 2001; Herrnstein & Ho 2003). 

In 2001, an upgrade of the 23 GHz system at the VLA enabled us to make the first 
interferometric map of the Galactic center in NH3(6,6) (z/ — 25.056025 GHz). With an 
energy above ground of 412 K, NH3(6,6) traces significantly warmer gas than our earlier 
NH3 data. An identical setup to our previous NH3 observations was used. After application 
of a Gaussian taper to the {u,v) data, the resulting synthesized beam is 16" x 14" (Her- 
rnstein & Ho 2002). A velocity integrated map of NH3(6,6) emission is shown in Figure 
Id. Contour levels are the same as in panels a-c, and correspond to steps of 4(j66, where 
(766 — 0.33 Jy Beam~^ km s~^. While some features such as the western streamer are also 
observed in lower NH3 rotation inversion transitions, the velocity integrated image is dom- 
inated by a cloud of hot molecular gas located less than 2 pc in projected distance from 
Sgr A*. This cloud is apparently located interior to the CND, and fills the molecular "hole" 
in the NH3(3,3) map (Herrnstein & Ho 2002). 



3. Derived Parameters 

3.1. Opacity and Intrinsic Line Width 

Hyperfine splitting of the NH3 rotation inversion transitions enables direct calculation 
of gas opacity (Ho & Towncs 1983). Each NH3 rotation inversion transition is split into 
five hyperfine components by the interaction between the electric quadrupole moment of the 
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nitrogen nucleus and the electric field of the electrons (see Townes & Schawlow (1975) and 
Ho & Townes (1983) for a detailed description of the NH3 rotation inversion transitions). 
Although spin-spin interactions further split each of these components, the resulting lines are 
closely spaced in frequency and arc only resolvable in cold, quiescent cores where line widths 
are less than 1-2 km s~^. Therefore, the spectral profile for a typical molecular cloud with 
a line width of a few kilometers per second will consist of a main line and two symmetric 
pairs of satellite hyperfine lines, which are located 10-30 km s~^ from the main line. At the 
Galactic center, however, molecular clouds have line widths in excess of 10 km s~^, and even 
the five electric quadrupole components are blended. The blending of these components 
produces a correlation between measured line widths of the gas and calculated main line 
opacities (McGary & Ho 2002). 

In McGary & Ho (2002), we presented a method to disentangle the opacity and intrinsic 
line width by using the observed line widths (Avobs) and main line flux densities of two 
rotation inversion transitions. In that paper, we used NH3(1,1) and (3,3) because they 
had the highest signal-to-noise (S/N) (see Figure 1). However, as a result of their different 
excitation requirements, the NH3(3,3) and (1,1) lines do not in fact trace the same volume of 
gas at the Galactic center (see §5.1). We have recalculated the NH3(1,1) main line opacity 
(Tni(l, 1)) and Avint using the NH3(1,1) and (2,2) data. With an energy difference of only 
42 K, NH3(1,1) and (2,2) are more likely to trace the same material. NH3(2,2) is significantly 
fainter than NH3(3,3), and the lower S/N hmits somewhat the number of pixels to which 
our method can be applied. 

In addition to changing the two rotation transitions considered in our method, we also 
use a Monte-Carlo simulation with 1000 trials to determine the significance of the derived 
parameter value at every pixel. Solutions are calculated following the algorithm of McGary 
& Ho (2002) for every pixel with a > 3(7 detection of both the (1,1) and (2,2) lines and 
an uncertainty in the observed line width of < 20 km s~^. Uncertainties in the measured 
flux densities are equal to the rms noise in the map, and uncertainties in the observed line 
widths are estimated from the fit of a Gaussian profile to the spectrum. Pixels included in 
the final map must successfully converge in at least 50% of the trials. The exact value of 
this cutoff does not appear to be important, and all pixels with solutions > 50% of the time 
show similar trends. 

For each pixel, the final value for the parameter is taken to be the mean value from the 
Monte-Carlo simulations. The mean values of Avobs and rm(l, 1) at each pixel are plotted in 
Figure 2a and c. No smoothing has been applied to the data, and regions of low S/N have a 
"pixelated" appearance. Uncertainties in each parameter are calculated from the standard 
deviation of the distribution of solutions at each pixel. Figures 2b and d plot the S/N of 
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Avint and Tin(l, 1), respectively, at each pixel. 

Table 1 lists the mean and standard deviation of derived parameters (including Avjnt 
and Tin(l, 1)) for the main molecular clouds in the central 10 pc. To avoid confusion with 
mean pixel values described above, we refer to the cloud-averaged value as the "characteris- 
tic value." The characteristic parameter value is calculated as the average of the parameter 
at all pixels where a value was determined from the Monte-Carlo simulations with > 3a 
significance. For the northern ridge and western streamer, some parameters had no deter- 
minations with > 3(7 significance. In these cases, all pixels were used in the calculation of 
the characteristic value and standard deviation. Inclusion of all pixels will tend to lower the 
characteristic value because smaller values will have a lower significance for a given absolute 
uncertainty. The characteristic value quoted in Table 1 is intended to provide a general value 
that can be easily compared between different clouds. However, it is important to note that 
many parameters show significant variations, such as gradients and local maxima/minima, 
within a single cloud. The amount of significant variation of each parameter within a cloud 
is characterized by the standard deviation, which is also listed in Table 1. 

Both the 50 km s'^ CMC (Aa = 120", AS = 60") and the central region of the southern 
streamer (30", —90") have NH3(1,1) main fine opacities significantly greater than one (see 
Figure 2a). The characteristic opacity that we derive for the southern streamer is very similar 
to the values measured by Coil & Ho (1999). The western streamer and northern ridge have 
significantly lower opacities with typical values of 0.3 — 0.4. Hatched regions in Figure 2b 
denote pixels for which only a lower limit for the opacity could be calculated. Although 
some of these values are the result of poor signal-to-noise, others may refiect real regions of 
high opacity. Many of these pixels have main line opacities near five, which corresponds to 
the maximum value of opacity considered in our solution method (McCary & Ho 2002). 

Intrinsic line widths in the central 10 pc tend to have values between 10 and 20 km s~^ 
(see Figure 2c). Significant increases in intrinsic line widths are observable near many cloud 
edges. However, the apparent random distribution of these features throughout the central 
10 pc makes it difficult to determine the exact way in which they formed. In general, these 
gradients may reflect dissipative effects at cloud edges possibly due to external pressures or 
tidal effects. For those pixels where only a lower limit on opacity can be calculated, the 
effect of blending on observed line width cannot be determined. The NH3(2,2) observed line 
width is less affected by blending of the hyperfine lines than Avobs(l, 1), and it is therefore 
used as the upper limit to the intrinsic line width for these pixels (Herrnstein 2003). 

Using the derived values of Tni(l, 1) and Avint along with the measured main line fiuxes, 
we can calculate the (2,2)-to-(l,l) rotational temperature, excitation temperature, and col- 
umn density, and ultimately estimate cloud masses. In the derivation of each parameter, we 
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perform a Monte-Carlo simulation to determine the mean value and associated uncertainties 
at each pixel. Pixels for which only a limiting value was calculated for Tin(l, 1) and Avint are 
not included in the following calculations. 



3.2. Rotational and Excitation Temperatures 

The ratio of main line emission from two NH3 rotation inversion transitions can be used 
to calculate the rotational temperature of molecular clouds at the Galactic center. Prom Ho 
& Townes (1983), the equation for rotational temperature calculated using NH3(1,1) and 
(2,2) is 

T ^ -41-5 K 

where ^Ta{.ii) is the ratio of emission from NH3(2,2) and (1,1) main hyperfine hues. (Be- 
cause of the similar frequencies of the NH3 transitions, AT^ depends linearly on flux density, 
and the ratio of main line flux densities can be used in place of the ratio of antenna temper- 
atures.) Using rm(l, 1) and the peak fluxes of the NH3(1,1) and (2,2) emission, we calculate 
NH3 (2,2)-to-(l,l) rotational temperatures. Mean values of Tr21 from the Monte-Carlo sim- 
ulations and the associated S/N estimates arc plotted in Pigurc 3a and b. In addition. Table 
1 lists the characteristic (2,2)-to-(l,l) rotational temperature for each cloud. 

(2,2)-to-(l,l) rotational temperatures in the central 10 pc are generally 25 K with un- 
certainties on the order of 10%. The relatively small uncertainties are the result of a weak 
dependence on rin(l, 1) for < rin(l, 1) < 1 and (2,2)-to-(l,l) line ratios less than one (Ho & 
Townes 1983). The derived rotational temperature depends primarily on the (2,2)-to-(l,l) 
line ratio, which is determined with a much higher S/N than Tm(l, 1). The coldest gas is 
found in the 50 km CMC (120", 60"), where the characteristic rotational temperature 
is 21 K. The western streamer (—80", —40" and —50", 40") has the warmest gas, with a 
characteristic (2,2)-to-(l,l) rotational temperature of 31 K. 

The measured rotational temperature is a lower limit on the true kinetic temperature 
of the gas (Martin et al. 1982). Table 1 hsts the kinetic temperatures calculated from 
{Tr2i) using the relations of Danby et al. (1988), which assume a H2 density of 10^ cm"^. 
(The degree to which Tr21 underestimates Tx decreases as the H2 density increases, until 
Tb21 ~ Tk for ~ 10^ cm^^ (Hiittcmcistcr et al. 1993a).) The resulting range of kinetic 
temperatures for clouds in the central 10 pc is roughly 20-50 K. It is important to note that 
NH3 (2,2)-to-(l,l) rotational temperatures arc not sensitive to gas with kinetic temperatures 
much greater than 50 K (see Danby et al. 1988). Therefore, the true range of temperatures at 
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the Galactic center is much larger than indicated by the (2,2)-to-(l,l) rotational temperatures 
in Figure 3a. This effect will be discussed in more detail in §5.1. 

Prom standard radiative transfer, the excitation temperature is related to the measured 
antenna temperature of the line emission, AT^, by 

ATa = v<^ {T,,-T,g){l-e-^) (2) 

where rj is the telescope efficiency (0.4 for the VLA at 23 GHz), $ is the filling factor, T^.^ 
is the excitation temperature of the gas, and T^g is the temperature of background radia- 
tion incident on the cloud (assumed to be the 2.7 K background of the Cosmic Microwave 
Background (CMB)). Solving for excitation temperature, we have 

Assuming the filling factor is equal to one (probably not a realistic assumption, see below), 
we have calculated the excitation temperature for gas in the central 10 pc. The resulting 
maps of mean excitation temperature and the S/N of these values determined from Monte- 
Carlo simulations are shown in Figure 3c and d. Overall, the excitation temperature has a 
smooth distribution (especially in regions with high S/N), with values ranging between 4 and 
8 K. The characteristic excitation temperatures that we derive are consistent with the peak 
values reported by Coil & Ho (1999, 2000). Our results are also consistent with excitation 
temperatures of ~ 5 K derived from NH3 observations within a few degrees of the Galactic 
center (Hiittemeister et al. 1993a). 

Throughout the central 10 pc, calculated excitation temperatures are much lower than 
the corresponding (2,2)-to-(l,l) rotational temperatures. There are two possible explana- 
tions for this discrepancy. First, low densities could lead to sub-thermal excitation due to 
lack of collisions. From detailed balance of emission and absorption, the equation for Tex is 
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^-hl//kbTex _ "crit 

n{H2) , , -• ' ^ 



where kf, is the Boltzmann constant and u = 1/ (e'*''/'^*-^'''"* — 1) (see e.g. Stutzki & Winnewisser 
1985). The critical density is defined as the ratio of the Einstein coefficient for spontaneous 
emission to the collision coefficient. Using the values for Aiq from Ho & Townes (1983) and 
collision coefficients from Danby et al. (1988), we calculate a critical density of 2 x 10^ cm~^ 
for every rotation inversion transition that we have observed. Figure 4 plots the excitation 
temperature as a function of H2 volume density for selected kinetic temperatures. Despite 
the low critical density, excitation temperatures do not increase significantly until n{H2) ~ 
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10^ cm~^, and the gas is not thermalized {T^x ~ Tk) until densities approach 10^ cm~^. The 

large differences between our measured excitation and kinetic temperatures imply an H2 
density of ^ 10"^ cm""^, which is significantly lower than densities measured for the GMCs 
and CND. Sub-thermal excitation of this sort would imply even larger kinetic temperatures 
for molecular gas in the central 10 pc. 

In §5.2, we will argue that at least the hottest gas components in the Galactic center 
must be nearly thermally excited with n{B.2) ^ 10^ cm~^. A more likely cause of the low 
measured excitation temperatures is small filling factors. If the clouds fill only a small fraction 
of the beam, then the observed excitation temperature will be lowered by the ratio of the 
cloud area to the area of the beam. If the molecular gas in the central 10 pc is thermalized, 
then the measured difference between T^.^ and Tr2i implies typical filling factors of ~ 0.2. 



The NH3(1,1) column density is related to rin(l, 1), Tex, and Avint by (Ho et al. 1977) 



By correcting for the fraction of molecules in the (1,1) state, which is characterized by the 
rotation temperature, A^n can be converted to the total NH3 column density, A^nhs (Townes 
& Schawlow 1975). Although emission from non-metastable states has been observed from 
the nearby Sgr B2 star-forming region (Hiittemeister et al. 1993b), we consider only the 
lowest 19 metastable states of NH3 in our calculation of the fractional population of the 
NH3(1,1) transition. The omission of non-metastable states from our calculation will at 
most cause an overestimation of the fractional population of a factor of two (Herrnstein 
2003). 

Total NH3 column densities determined from our data are plotted in Figure 5 and 
characteristic values for specific clouds are listed in Table 1. Characteristic column densities 
for the main features range from 2 — 9 x 10^^ cm^. However, the largest column densities 
(A^NHs > 10^^ cm^) are found primarily in the features associated with the GMCs. The central 
core of the southern streamer (30", —80") has a NH3 column density of 1.4(ih0.3) x 10^^ cm~^. 
The 50 km s~^ GMC also has large NH3 column densities, with highest values found at 
the northeastern edge of our mosaic. This GMC extends more than an arcminute past the 
northeastern edge of our map (see Figure 9 in McGary et al. (2001)), and emission at (—130", 
80") is associated with dense regions of the cloud where column densities should be highest. 



3.3. 



Column Densities and Masses 
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At the Galactic center, the total mass of a cloud is given by 

M = 2.4 M, ( ,) (^) (■^) " . (6) 

where A is the area of the cloud in arcscc^ and X(NH3) is the abundance of NH3 relative to 
H2. At present, X(NH3) is not well-known for Galactic center molecular clouds. Abundances 
measured in the Perseus and -3 km s^^ spiral arms and in dense, star-less cores suggest 
A:(NH3)= 10-^ - 10-^ (Benson & Myers 1983; Batrla et al. 1984; Serabyn & Guesten 1986). 
However, there is much evidence that NH3 abundances are enhanced in warm clouds, most 
likely as the result of evaporation of NH3 molecules off of dust grains (Pauls et al. 1983; 
Walmsley et al. 1987). NH3 abundances range from 10"^ to 10~^ in warm molecular clouds, 
with the highest measured abundances coming from the Sgr B2 complex, one of the largest 
star formation regions in the Galaxy (Hiittemeister et al. 1993b, and references therein). The 
precise relationship between X(NH3) and is not straightforward. Assuming that high 
temperatures at the Galactic center will produce enhanced abundances of NH3, we adopt 
X(NH3)=10~^ in our calculations. Resulting estimates of masses are presented in Table 1. 
The assumption of an enhanced NH3 abundance should produce smaller cloud masses than 
those calculated by Coil & Ho (1999, 2000), who use ^(NHg)^ 10"^ 

A comparison of our derived mass for the edge of the 50 km s~^ GMC to published 
values in the literature can be used as a consistency check for our choice of X(NH3). Total 
mass estimates for the 50 km s~^ GMC range between 6 x 10^ and 4 x 10^ Mq (Giisten et al. 
1981; Armstrong & Barrett 1985; Zylka et al. 1990). The virial relation between the mass 
of a cloud and its size and line width (My 2 x 10^ M© (^r^j ( 1 ^ density 

profile of 1/r, MacLaren et al. (1988)) can also be used to calculate the minimum mass of 
a gravitationally bound cloud. Assuming a size of 5 pc and a line width of 17 km s"^ (see 
Table 1), the virial mass of the 50 km s^^ GMC is ~ 3 x 10^ Mq. As previously mentioned, 
roughly 3/4 of the 50 km s~^ GMC, including the densest condensations, lies outside our 
VLA mosaic. Our derived mass of 5 x 10^ Mq for the fraction of this cloud within our 
mosaic is therefore consistent with what we expect, and we conclude that X(NH3)= 10"^ is 
a reasonable assumption for molecular clouds near Sgr A*. 

The southern streamer is the second-most massive feature in the central 10 pc. We 
estimate a total mass of 8 x 10*^ Mq for the southern streamer. Using X(NH3)= 10^^, Coil 
& Ho (1999) calculate a mass of 6 x 10"^ Mq for gas from the "central cloud" of the southern 
streamer, which is roughly coincident with the core of bright NH3 emission centered at (-1-30", 
—80") in our maps. If the Coil & Ho (1999) result is re-calculated using our assumed value 
of X(NH3) then the associated mass would be 6 x 10^ Mq, more than an order of magnitude 
smaller than our value. The large difference between these results reflects different definitions 
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for the size and extent of the southern streamer. First, Coil & Ho (1999) assume a size of 
only 300 arcsec^ for the central cloud, while it appears to be closer to 2000 arcsec^ in our 
velocity integrated maps. In addition, the emission associated with the southern streamer at 
(0", —125") is associated with the northern tip of the "southern cloud" in Coil & Ho (1999) 
and not included in the above mass estimate. It should also be noted, that the measured NH3 
column for the central cloud is consistent between the two papers (for the same assumption 
of X(NH3)). Considering the different assumptions outlined above, the reported masses for 
the southern streamer are consistent. However, we believe that our result is a more complete 
measure of the mass of the southern streamer from 50" to 150" south of Sgr A*. 

We estimate masses of 1 x 10"^ Mq and 2 x 10^ Mq for SEl and SE2, respectively. 
Together, these clouds have a mass of a few x 10^ Mq, and constitute the majority of the 
dense "molecular ridge" that connects the 20 km s"^ CMC in the south to the 50 km s"^ 
CMC in the northeast (Coil & Ho 2000). Both the western streamer and northern ridge have 
masses roughly an order of magnitude smaller than the southern streamer and two orders of 
magnitude smaller than estimated total masses of the 20 and 50 km s~^ GMCs (Zylka et al. 
1990). As we argue in the following section, we believe these smaller clouds have been swept 
up by the expansion of Sgr A East. 

4. The Impact of Sgr A East on Nearby MoleculEir Clouds 

Maps of radio continuum emission from the Galactic center show a large shell of syn- 
chrotron emission with a radius of ~ 1.6' (4 pc) centered roughly 1' to the east of Sgr A* 
(Pedlar et al. 1989). This shell is called Sgr A East and is thought to lie within a few parsecs 
of the Galactic center (see §7). To the northeast, the edge of Sgr A East lies precisely along 
the western edge of the 50 km s~^ GMC, and it is generally agreed that Sgr A East is in- 
teracting with the GMC to some extent. Previous work by Mezger et al. (1989) interpreted 
this interaction as indicating that much of the material in the 50 km s~^ GMC was cleared 
out of the central parsecs by the expansion of Sgr A East. The resulting estimate of the 
energy of Sgr A East is ^4x10^^ ergs, more than an order of magnitude larger than typical 
supernova remnants (Mezger et al. 1989). This result has sparked many discussions on the 
possible origin of Sgr A East, including multiple correlated supernovac, expansion within 
a bubble (Mezger et al. 1989), or tidal disruption of a star by the supermassive black hole 
(Khokhlov & Meha 1996). 

The assumption that a significant amount of material in the 50 km s~^ GMC was swept 
up by the expansion of Sgr A East was based primarily on the apparent concave morphology 
of the GMC, with the northern part of the cloud wrapping around the top of the Sgr A East 
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shell (Mezger et al. 1989). However, in McGary et al. (2001), we showed that this "northern 
ridge" is kinematically distinct from the 50 km s""*^ GMC, with a characteristic velocity of 
-10 km s~^. Using the derived cloud parameters presented in the previous sections, we 
no longer need to rely on morphology as the primary indication of an interaction. In the 
following paragraphs, we revisit the question of the impact of Sgr A East on molecular gas 
near Sgr A*, and we calculate a new estimate for the energy associated with the progenitor 
explosion. 

The 50 km s^^ GMC, northern ridge, and western streamer lie along the edge of 
Sgr A East (McGary et al. 2001). However, the physical conditions differ significantly be- 
tween these clouds. The molecular gas in the 50 km s^^ GMC does not appear to be strongly 
affected by the impact of Sgr A East. The characteristic (2,2)-to-(l,l) rotational tempera- 
ture of the 50 km s~^ GMC {{Tr2i) — 21 K) is the lowest value of any feature in the central 
10 pc. If the 50 km GMC is moving outwards with the expansion of Sgr A East, then 
we should expect to see a velocity gradient along the cloud. However, no velocity gradient 
is observed in this cloud (McGary et al. 2001). 

The northern ridge and western streamer are more than a magnitude less massive than 
the 50 km s~^ GMC, and thus might be expected to be more affected by the impact of 
Sgr A East. The northern ridge is a linear feature, roughly 90" in length, that lies along 
the northern edge of Sgr A East. Compared to the 50 km s^^ GMC, the northern ridge has 
a similar characteristic intrinsic line width ((Avint) = 17 km s^^) and higher characteristic 
(2,2)-to-(l,l) rotational temperature ((Tr2i) = 28 K). There is no velocity gradient larger 
than the size of the spectral resolution in our data (~ 9.8 km s^^), but this result is not sur- 
prising because the cloud velocity of -10 km s~^ is consistent with the bulk of any expansion 
occurring perpendicular to the LOS. 

The strongest evidence for an interaction between Sgr A East and molecular gas comes 
from the western streamer. The western streamer also has a long, filamentary structure, 
and its curvature matches closely the western edge of Sgr A East. In McGary et al. (2001) 
we reported a velocity gradient of 1 km s~^ arcsec"^ (25 km pc~^) covering 150" (6 pc) 
along the western streamer. At (60", —80") emission is at —70 km s~^, but the velocity 
smoothly shifts towards the red until it reaches ~ 90 km s~^ at (—70", 70") (McGary et al. 
2001). A velocity gradient of this type is consistent with a ridge of gas highly inclined to 
the line-of-sight that is being pushed outwards by expansion of Sgr A East. In this model, 
the southern part of the streamer would be located on the front side of the shell, thereby 
producing the observed blue-shifted velocities. Molecular gas in the western streamer has 
the largest characteristic intrinsic fine width (20 km s"^) and highest temperatures in the 
central 10 pc. The characteristic (2,2)-to-(l,l) rotational temperature of 31 K corresponds 
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to a kinetic temperature near 46 K. The western streamer is also the only one of these three 
clouds with strong NH3(6,6) emission, and it is hkely that some fraction of the gas has kinetic 
temperatures ^ 100 K. 

Based on the above comparison, it seems unlikely that a significant amount of the 50 
km s~^ GMC has been moved by Sgr A East. However, there is convincing evidence that the 
expansion of Sgr A East swept up the material in the western streamer. Assuming that only 
the western streamer was cleared out of the central parsecs by Sgr A East, we calculate a 
new estimate of the age and energy of the progenitor explosion in the following paragraphs. 

A supernova remnant (SNR) is reckoned to have entered the isothermal or "snow-plow" 
phase when the time scale for radiative cooling becomes less than the dynamical time of the 
SNR. At this point, the shell expands with constant radial momentum, and a thin, dense 
layer of gas is formed that moves outwards with the shock front. For a SNR in the snow- 
plow phase, the energy of the progenitor supernova, Esn, is related to the current energy 
associated with the shell, Egh, by (ShuU 1980) 



Esh J \R. 
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where Rsu is the radius of the shell and Rgg is the radius at the time of shell generation, tsg- 
The shell radius at time, t, is given by (ShuU 1980) 
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where no is the mean initial H density in cm~^. At t = tsg-, Equation 8 defines the shell 
generation radius, Rsg. Substituting the equation for Rgg into Equation 7, and solving for 
Esn we find 
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Esn = 2.0 x 10^° ergs ( Vj^ V ( \ ^ _ 

^ V 10^ ergs/ \lpcj Vl00cm-3y ^' 

If the 4 X 10^ Mq of material in the western streamer was swept out of a region with 
a radius of 4 pc, then the implied mean initial density is Uq ~ 600 cm~^, which is typical 
of molecular clouds. The kinetic energy associated with the western streamer (assuming a 
velocity of 100 km s"^ (McGary et al. 2001)) is E^s = 4 x 10^° ergs. Taking as the shell 
energy, the derived energy for the progenitor supernova explosion is Esn = 2 x 10^-*^ ergs. 
This estimate of Esn is a lower limit on the energy of the progenitor explosion because the 
western streamer covers only a fraction of the total shell. The shell energy (and therefore 
Esn) will scale as 47r/fi^s, where corresponds to the solid angle initially covered by 
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the material in the western streamer. This sohd angle likely differs significantly from the 
solid angle that it currently covers. For example, as the shock front expanded, it may have 
wrapped around the densest regions compressing them into narrow filaments and thereby 
decreasing the solid angle of the material in the western streamer. Based on observations 
of ionized material surrounding Sgr A East, the total kinetic energy of the shell has been 
estimated to be ~ 2 x 10^^ ergs (Mezger et al. 1989; Genzel et al. 1990). This kinetic energy 
would imply that Jl^s/47r ~ 0.15. The mean initial density calculated from the western 
streamer is also a lower limit. However, the observed distribution of dense molecular gas 
in the region indicates that Sgr A East probably expanded into a clumpy and highly non- 
homogeneous environment. The filaments of dense gas detected in NH3 emission therefore 
likely account for the majority of the gas mass surrounding the remnant, and the estimate 
of rio from the western streamer is probably accurate to within a factor of a few. Taking 
Esh = 2 X 10^^ ergs and assuming rio = 1000 cm^"^, we calculate an upper limit on the energy 
of the progenitor supernova of at = 9 x 10^^ ergs. 

Based on the range of energies for the progenitor explosion calculated from our molecular 
data, we conclude that Sgr A East is most likely the result of a single supernova that occurred 
near the Galactic center. The expansion of Sgr A East into molecular gas produced the 
western streamer and possibly the northern ridge, but it has had little effect on the majority 
of the material in the 50 km s~^ GMC. In fact, this GMC is likely significantly slowing the 
expansion of the supernova remnant to the northeast of Sgr A*. 

The case for a single supernova explosion as the progenitor of Sgr A East has been 
gaining momentum in the past few years. In particular, the small gas mass and thermal 
energy of 10^^ ergs inferred from Chandra observations of hot gas associated with Sgr A 
East is consistent with the ejecta of a single supernova (Maeda et al. 2002). Based on the 
centrally concentrated X-ray emission, Maeda et al. (2002) classify Sgr A East as a "mixed 
morphology" (MM) SNR produced by a Type II supernova. This classification is supported 
by the detection of multiple 1720 MHz OH mascrs, which are often found in MM SNRs, 
along the edge of the shell (Yusef-Zadeh et al. 1999; Green et al. 1997). 

Using the estimated energy of the progenitor supernova, we can calculate the age of Sgr 
A East. From ShuU (1980), the time of shell generation is given by 



Substituting the equation for tsg into Equation 8, we determine an age of (1 — 2) x 10^ yr 
for Sgr A East. This age is similar to the age of 7500 yr derived by Mezger et al. (1989) 
for a model in which the supernova occurred within a low-density bubble inside the dense 
50 km s~^ GMC. However, the lack of evidence for a strong impact on the 50 km s~^ GMC 
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leads us to favor a scenario in which the supernova occurred outside the GMC. Our derived 

age of ~ 10^ yr agrees well with other independent measurements of the age of Sgr A East. 
Based on the expected differential shearing of a shell near the Galactic center, Uchida et al. 
(1998) independently derive an age of a few 10^ year for Sgr A East. An age of 10^ yr is also 
consistent with the X-ray properties of the shell (Maeda et al. 2002). 

As noted elsewhere (see e.g. Serabyn et al. 1992; Morris & Scrabyn 1996; Maeda et al. 
2002) the inferred age of the supernova is 1-2 orders of magnitude smaller than the age of 
the Hll regions that lie along the western edge of the 50 km s~^ GMC. Therefore, these Hll 
regions could not have been triggered by the impact of the Sgr A East shock front on the 
molecular cloud. However, the alignment of these star formation regions along the edge of 
the 50 km s"^ GMC is noteworthy, and may be the result of previous interactions between 
this cloud and expanding shock fronts that originated near the nucleus. 



5. The Temperature Structure of the Central 10 pc 



5.1. A Two- Temperature Gas 

In general, a rotational temperature can be calculated from any two rotation inversion 
transitions of NH3 using the equation 
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where A£'„/ is the energy difference between the upper {J,K — u) and lower {J,K — I) 
rotation inversion states and Cui = i''u!°|2fc relates the statistical weights {gj), magnetic 



dipole moment (|/Xio|i), fraction of emission in the main line (/cj) of the two transitions. 
This equation assumes equal excitation temperatures, filling factors, and telescope efficiencies 
for the two transitions. Due to slow mixing rates (~ 10~^ yr~^) between the "ortho" {K = 
3n) and "para" {K ^ 3n) spin states of NH3, rotational temperatures are generally calculated 
using two transitions from the same spin state (Ho & Townes 1983). Therefore, if the gas is 
in LTE, then we can calculate a (6,6)-to-(3,3) rotational temperature from our VLA data. 

NH3 (6,6)-to-(3,3) line ratios in the central 10 pc are higher than would be expected 
based on NH3 (2,2)-to-(l,l) rotational temperatures (Herrnstein & Ho 2002). Within 2 pc of 
Sgr A*, the (6,6)-to-(3,3) hue ratios exceed theoretical hmits, preventing calculation of Tr63 
(Herrnstein & Ho 2002; Herrnstein 2003). These "un-physical" line ratios will be explained 
in §5.2. Elsewhere in the map, faint NH3(6,6) emission makes calculation of Trgs difficult 
(see Figure 1). In addition, the NH3(3,3) main line opacity cannot be reliably calculated 
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from Tin(l, 1) because these two transitions likely trace gas at different temperatures (see 
below). Although rin(3, 3) is directly related to the ratio of the main and outer-most satellite 
hyperfine lines, the relative population of the satellite lines compared to the main line is only 
0.008 (compared to 0.22 for NH3(1,1)), (Townes & Schawlow 1975)), and emission from the 
satellite hyperfine lines is generally below the rms noise level of our map. 

Due to the above limitations, T/jg3 can only be directly calculated for the core of the 
southern streamer and SEl. At (33", —75"), the southern streamer has a measured NH3 
(6,6)-to-(3,3) line ratio of 0.16 ± 0.02 and a NH3(3,3) main line opacity of Tiii(3, 3) — 2.5^*^0^3. 
The infcrrcd(6,6)-to-(3,3) rotational temperature of Trss = 80±5 K corresponds to a kinetic 
temperature of lOOl^g K for gas with a H2 density of 10^ cm~^ (Danby et al. 1988). In SEl 
(73", —87"), the NH3(3,3) main line opacity is 3.6l|^;5. The (6,6)-to-(3,3) main line ratio of 
0.27lo;5 corresponds to Trss = 88+^4 K and a kinetic temperature of IIO+20 

Lower limits on the kinetic temperature can be calculated in SE2 and the 50 km s~^ 
GMC. In SE2 (108",-60"), the NH3(6,6)-to-(3,3) line ratio is 0.39l[J:|Ji and an upper limit of 
1.5 is calculated for Tin(3, 3). These values result in lower hmits for the rotational and kinetic 
temperatures of Trgs ~ 120 K and Tk ^ 190 K. For the 50 km GMC, NH3(6,6) is only 
detected near (112", 48"). (Although significant, the velocity integrated NH3(6,6) emission at 
this position does not exceed 1.32 Jy beam~^ km s~^ 4o"66) and thus is not seen in Figure 
Id.) At this position, the (6,6)-to-(3,3) line ratio is 0.16±0.05 and rm(3, 3) ^ 2. The inferred 
lower limits on the derived temperatures are T^qs ^ 75 K and ^ 90 K. Due to low S/N 
and large line widths, the main line opacity cannot be estimated for the western streamer. 
The typical line ratio of 0.7 in this cloud implies an upper limit of Trgs ~ 230 K (Herrnstein 
& Ho 2002). This upper hmit on the (6,6)-to-(3,3) rotational temperature provides little 
constraint on the kinetic temperature of the gas (Danby et al. 1988). 

The kinetic temperatures derived above are significantly higher than those derived from 
NH3 (2,2)-to-(l,l) line ratios. This discrepancy is indicative of a gas composed of a mixture 
of separate cool and hot components. If multiple temperature components are contained 
within the telescope beam, then emission from different rotation inversion transitions will not 
trace the same material, and the resultant measurement of kinetic temperature will depend 
strongly on the choice of rotation inversion transitions used in the calculation. Because the 
fractional population of the (J,K) state is a function of kinetic temperature, each NH3(J,K) 
transition is most sensitive to gas at a particular temperature corresponding to the peak 
fractional population of that transition. The peak fractional population of NH3(1,1) and (2,2) 
occurs around 20-30 K (see Herrnstein 2003), making these transitions excellent tracers of 
cool molecular material. At kinetic temperatures greater than 50 K, however, Tr2i becomes 
insensitive to changes in Tk- All kinetic temperatures ^ 80 K have (2,2)-to-(l,l) rotational 
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temperatures between 40 and 60 K (Walmsley & Ungerechts 1983; Ho & Townes 1983; Danby 
et al. 1988). Higher transitions such as NH3(3,3) and (6,6) have a small fractional population 
in cool clouds. (For example, only ~ 1 % of the NH3 molecules are in the NH3(6,6) state 
for gas at 50 K.) The fractional population rapidly increases towards higher temperatures, 
making it easier to detect the emission. In addition, Tr63 is much more sensitive to variations 
in kinetic temperatures between 50 and a few hundred Kelvin. While an increase in kinetic 
temperature from 100 to 200 K changes Tr2i by only a few degrees, Tr63 increases from 80 
to 125 K. 

Extensive singlc-dish observations of NH3 rotation inversion transitions by Hiittemeister 
et al. (1993a) show that molecular clouds in the central 4° x 1° of the Galaxy are best 
characterized by a two-temperature gas distribution on the scale of 3 pc (corresponding to 
the size of the NRAO 40-m primary beam). In general, 75% of the total column of NH3 
comes from a cool gas at 25 K, while the remaining 25% of the material is contained in 
a hot component with a temperature of roughly 200 K (Hiittemeister et al. 1993a). The 
uncertainty in these temperatures is roughly 30%, but the cold temperature matches well 
the observed dust temperature in the region. 

Based on the temperatures calculated in this paper, it appears that a similar cloud 

structure is valid in the central 10 pc on a spatial scales of ~ 0.5 pc. In Figure 6a, we 
plot the inferred kinetic temperature as a function of rotational temperature for different 
combinations of hot and cold gas. The hot gas is assumed to be at 200 K based on the results 
of Hiittemeister et al. (1993a). The cold component is set at 15 K in order to be consistent 
with the kinetic temperatures inferred from our Tr2i measurements. Models for the derived 
rotational temperatures are plotted for cases in which the hot component comprises 0, 25, 
50, 75, and 100% of the gas. 

Mean kinetic temperatures derived from our data arc overlaid in Figure 6a. The mean 
kinetic temperature derived from Tr2i is 26 ± 2 K. The cold component therefore must have 
a kinetic temperature lower than ~ 25 K. The two measurements of Trq^ from our data 
have a weighted average of 101 ± 9 K. Although the low measurements from 7^21 imply a 
lower temperature for the cold component than derived by Hiittemeister et al. (1993a), the 
data are consistent with their conclusions that roughly 25% of the gas is contained in a hot 
component. 

In order to study the temperature structure in the central parsecs in more detail, we 
have compiled pubhshed NH3 rotational temperatures from a 7.5' x 11' (17.5 x 25.5 pc) 

region surrounding Sgr A* that includes both the 20 and 50 km s"-*^ GMC. A summary of 
the references for the measurements, including the telescope and corresponding beam size, 
is given in Table 2, and the specific values for the rotational temperatures are listed in Table 
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3. Only one pointing from the large-scale survey by Hiittemeister et al. (1993a) is contained 
within this region. 

Figure 7, overlays the position of each observation on a 1.2 mm continuum image showing 
bright dust emission from the GMCs and free-free emission from ionized arcs of Sgr A 
West (Zylka et al. 1998). For single-dish measurements (Figure 7a), the size of the circle 
corresponds to the FWHM telescope beam. Single-dish observations have focused primarily 
on the large GMCs, and pointings tend to be arranged parallel to the Galactic plane. The 
FWHM of the primary beam is plotted for each interferometric observation in Figure 7b. 
Both Coil & Ho (1999, 2000) and the data presented in this paper include a pointing centered 
on Sgr A*. For the interferometric data, independent temperature measurements may be 
made once per synthesized beam. However, Okumura et al. (1989) and Coil & Ho (1999, 
2000) present calculated temperatures at selected positions only. We plot the synthesized 
beam for the NMA at the positions of rotational temperature measurements in Okumura 
et al. (1989). In Coil & Ho (1999, 2000), characteristic rotational temperatures are given for 
the main condensations and we report their estimated positions in Table 3. 

Using the relation of Danby et al. (1988), kinetic temperatures have been calculated 
for all measurements of Tr2i, 7ii42, 7ii54, and Tr63- (2r4i is not included in Danby et al. 

(1988), and it is therefore not converted to a kinetic temperature.) In Figure 8, we plot 
kinetic temperatures as a function of distance from Sgr A*. The uncertainty in position due 
to the size of the beam (see Table 2) is not plotted. We also calculate a radial profile for 
the kinetic temperatures derived from our Tr21 measurements. The best fit radial profile 
[Tk = 31(±2) — 0.06(±0.02)r] is overlaid as the grey line in Figure 8a. The surrounding 
light grey box shows the approximate range of derived kinetic temperatures from pixels with 
a > 3(7 determination of Tji2i. 

In general, kinetic temperatures in Figures 8a, b, and d are consistent with a relatively 
constant kinetic temperature within 500" of Sgr A*. Kinetic temperatures derived from our 
Tr2i measurements show a slight increase as gas nears Sgr A*. Unfortunately, there are no 
measurements of Tr63) which would be more sensitive to large increases in the kinetic tem- 
perature, within 100" of the nucleus. Kinetic temperatures derived from 7^63 for distances 
between 100" and 150" are consistent with measured kinetic temperatures at larger radii. 
Any significant heating of the molecular gas is therefore confined to a region less than ~ 4 pc 
from Sgr A *. 

In Figure 6b, the average kinetic temperatures calculated from Tji2i, 3ij42, and Tr63 
measurements in the central 20 pc are plotted as filled circles overlaid on the same models 
for a two-temperature gas as Figure 6a. These temperatures are calculated as the average of 
the measurements in Figure 8 weighted by the associated uncertainties (see Table 3), and are 
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equal to 26 ±6 K, 55 ±5 K, and 100 ±8 K, for kinetic temperatures derived from Tr2i, 7ii42, 
and 7ij63, respectively. (Values of Tr2i reported by Glisten et al. (1981) differ significantly 
from other published rotational temperatures (see Table 3) and imply extraordinarily high 
kinetic temperatures. Omission of these data from our analysis results in (Tr2i) = 24 it 2 K 
and does not affect the conclusions below.) The mean kinetic temperatures are consistent 
with the kinetic temperatures measured from our VLA data alone. Based on the consistency 
between our VLA measurements and these previously published results, we conclude that 
molecular clouds within 10-20 pc of Sgr A* have a two-temperature structure on size scales 
of ^ 0.5 pc in which roughly one quarter of the gas is contained in a hot component. 

Only one measurement of Tr54 has been made near the 20 or 50 km s^^ GMC (Hiittenieister 
et al. 1993a). The inferred kinetic temperature from this measurement is marked by a tri- 
angle in Figure 6b. The corresponding kinetic temperature derived from Tr2i at the same 
position is also plotted. Both kinetic temperature measurements are significantly larger than 
expected for our two-temperature model. However, there is a large scatter in Tr2i measure- 
ments in this region, and it is not surprising that this single measurement set does not follow 
the same trend as the weighted average result. 



5.2. Interpreting the Highest Line Ratios at the Galactic Center 
Rotational temperatures calculated using Equation 11 will give a reasonable value only 

if 

< Z^ln [1 - R^^^i (1 - e-^-O] < 1 ' (12) 

where Rm,ui = \"'o^ observed main line ratio of the upper and lower transitions. 

This constraint translates to an upper limit on the main line ratio of 

Rm,ui < — --. 3;-- — ) (13) 

1 — e 

which has a maximum value of 1/Cui when rm,i <^ 1- Thus 1/C„/ is a firm upper limit on the 
observed line ratio, regardless of the opacity or kinetic temperature of the gas. We refer to 
hne ratios that exceed this limit as "un-physical" because they are not theoretically possible 
for a single cloud in LTE. 

It is important to note that the two-temperature model described in the previous sec- 
tion cannot account for un-physical line ratios. In that model, gas components at two 
different temperatures, but each in LTE, are contained within the telescope beam. Because 
the observed line ratio is simply the weighted average of the line ratios from the two gas 
components, it will never exceed the line ratio associated with the hot cloud component. 
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Un-physical (6,6)-to-(3,3) line ratios (where Rm,63 > 2.3) are found throughout the 

central 2 pc. Figure 9 shows the positions of 12 spectra overlaid on velocity integrated 
NH3(6,6) emission. Figure 10 plots the corresponding NH3(3,3) (grey) and NH3(6,6) (black) 
spectra. The spectra have been Hanning smoothed, and have an rms noise in one channel 
of 7.6 and 9.5 mJy for NH3(3,3) and (6,6), respectively. Spectrum A is coincident with Sgr 
A*. In this section, we focus on the spectra associated with un-physical line ratios. Gas 
kinematics will be discussed in the following section. 

Un-physical NH3 (6,6)-to-(3,3) line ratios occur near -20 km s^^ in spectra A and B, and 
-|-70 km s^^ in spectra D and E. These line ratios do not appear to result from an enhancement 
in NH3(6,6) emission. Instead, they occur at velocities where there is no NH3(3,3) emission. 
This complete lack of NH3(3,3) emission implies that un-physical line ratios are likely the 
result of an absorption process along the LOS. The most straightforward scenario for the 
origin of un-physical line ratios is one in which emission from a hot cloud with a high (but 
physical) (6-6)-to-(3,3) line ratio, passes through cold molecular material along the line of 
sight. The lower-energy NH3 (3,3) emission will be significantly absorbed by the cold cloud, 
while NH3(6,6) will be relatively unaffected. This model is drawn schematically in Figure 
11 and explored in more detail in the following paragraphs. 

For simplicity, we assume that there is no radiative coupling between the two clouds in 
our model. From basic radiative transfer, the brightness temperature of the first cloud is 
expressed as 

n,i = Te,,i(l - + Ttge-^' . (14) 

For the cold cloud we assume that the background temperature is equal to the temperature 
of the CMB. If Cloud 2 is near Cloud 1, then emission from Cloud 1 will dominate the 
background, and the brightness temperature of the second cloud is 

n,2 = 7;,,2(1 - e-^^) + Tfe.ie-^^ . (15) 



The antenna temperatures "on" and "off" source are related to the brightness temper- 
ature of Cloud 2 by 

T4,on = [<I>T,,2 + (1 - $)Tcmb] (16) 
Ta,oS = rjTcMB , (17) 

where rj is the telescope efficiency and $ is the fiUing factor of the source in the telescope 
beam. The difference between these two antenna temperatures will correspond to the mea- 
sm^cd antenna temperature due to the source (AT^). In practice, AT4 must be corrected 
for the fraction of emission in the main hyperfine component {k). For NH3(3,3) and (6,6), 
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however, the population of satellite hyperfine lines relative to the main component is very 
small and can be ignored {kg/k-^ = 1.08). Assuming equal efficiencies and filling factors for 
NH3(3,3) and (6,6), the observed main line ratio is therefore 



Ar^(6,6) ^ Tb,2(6, 6) - TcMB 



(18) 



In order to produce un-physical line ratios, Cloud 2 must have r2(3,3) ^ 1 and 
T2{Q, 6) <^ 1. This can be accomplished most efficiently for a cold cloud with Tc ~ 30 K. Ki- 
netic temperatures < 30 K have only minimal populations in both transitions, and although 
the relative population of NH3(3,3) compared to NH3(6,6) is higher at temperatures lower 
than 30 K, the small population of these transitions requires a large H2 column density. 
Clouds with temperatures greater than 40 K will begin to have a significant population of 
NH3(6,6) and both lines will be affected by absorption in the foreground cloud. 

Assuming r2(3, 3) 3> 1, the brightness temperature of the second cloud will be approx- 
imately equal to rex,2(3,3). For NH3(6,6), however, T2(6, 6) <^ 1 implies that T6_2(6,6) 
21),i(6, 6). Therefore, the observed line ratio can be approximately written as 

P _ 7b,i(6,6) -TcMB .^px 

-t ex,2l<J, <JJ — -t CMB 

Un-physical line ratios with Rm,63 ^ 1 imply that T6i(6,6) ^ Tex-2(3,3). Assuming that 
the NH3(6,6) excitation temperature in Cloud 1 exceeds Tcmb, then Equation 16 implies 
that Tf,^i(6,6) < Tea;,i(6,6). Therefore, un-physical line ratios will generally imply that 
T'ex,i(6,6) ^ Tej;,2(3, 3). As seen in §3.2, such large differences in excitation temperatures 
imply that the H2 density must exceed 10^ cm~^ in these clouds (see Figure 4). 

It should be noted that un-physical (6,6)-to-(3,3) fine ratios were not reported in the 
survey of Galactic center molecular clouds by Armstrong & Barrett (1985). However, these 
data do not include any molecular gas within 2 pc of Sgr A* (see Figure 8d). In addition, 
the single-dish measurements have a resolution of only 84", and thus sample the average 
temperature structure on scales of ~ 3 pc. Clouds with high or un-physical line ratios are 
likely to be small and therefore fill only a fraction of the telescope beam. Cool gas, however, 
is prevalent throughout the entire region. The inclusion of significant amounts of cool gas 
within the telescope beam will lower the observed line ratios. Because our VLA data have 
a resolution of 16" x 14", they are more sensitive to the temperature structure of molecular 
clouds on small scales. However, evidence for beam dilution of high line ratios is still detected 
at velocities near 20 km s~^. Although some (6,6)-to-(3,3) line ratios greater than one occur 
near 20 km s~^ (see Spectra C, D, and E), no line ratios exceed ~ 2. It is likely that any un- 
physical line ratios at these velocities are diluted by emission from cool gas in the 20 km s~^ 
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GMC. If beam dilution greatly affects observed line ratios, then the cloud that produces the 
un-physical line ratios must be significantly smaller than the synthesized beam of our VLA 
data. 

It is also worth noting that the maximum allowed line ratio can exceed 1/Cui if the 
cloud is not in LTE. Equation 11 assumes equal telescope efficiencies, fiUing factors, and 
excitation temperatures for the two transitions. If the values of these parameters for the 
upper transition exceed those for the lower transition, then the corresponding limit on the 
line ratio will exceed 1/C„/. However, because the energy separation of the inversion doublets 
is very small, it seems difficult to imagine a scenario in which the excitation temperatures 
of the two transitions differ significantly. We prefer the absorption model outlined above 
because it provides a simple explanation for un-physical line ratios that is consistent with 
observed spectral line profiles without having to invoke non-LTE conditions. 

In summary, we believe that the un-physical line ratios observed to the east and south- 
east of Sgr A* are the result of absorption of NH3(3,3) emission from a hot cloud near Sgr A* 
by cool material along the LOS. Although the current resolution of these data prevents a 
determination of the exact parameters of this cloud, the un-physical line ratios require that 
at least the hot material must be contained in a dense component with ~ 10^ cm~^. 

6. Sepeirating Multiple Gas Components Near Sgr A* 

Spectra taken from the central 2 pc of the Galaxy are very comphcated, with strong 
evidence for multiple clouds projected along the LOS (see spectra A-F in Figure 10). The 
(6,6)-to-(3,3) line ratios can be used to separate emission from the absorbed hot cloud from 
emission associated with cool, un-absorbed gas. Direct comparison of line ratios can be 
difficult because many features have infinite NH3 (6,6)-to-(3,3) line ratios. We therefore 
separate these gas components by using the difference between the NIl3(3,3) and (6,6) emis- 
sion. Gas with 5',^(6,6) > S';^(3,3) is termed "high line ratio" (HLR) gas, while gas with 
Sy{Q^ 6) < Sy{?), 3) is called "low line ratio" (LLR) gas. This designation works well because 
gas with 5*1,(6, 6) = S^if^i 3) (or .^^,63 = 1) corresponds to a rotational temperature of 340 K, 
or a kinetic temperature of thousands of Kelvin. Because most NH3 molecules will be de- 
stroyed at these temperatures, (6,6)-to-(3,3) line ratios greater than one most likely result 
from absorption of NH3(3,3) emission from a hot cloud by cool intervening material. 

The difference in NH3(3,3) and (6,6) flux is calculated at every pixel {x,y,v) where the 
emission exceeds three sigma in at least one of the tracers. Figure 12 plots velocity integrated 
NH3(6,6) emission associated with each gas component. The left-hand panel of Figure 12 
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plots velocity integrated map of NH3(6,6) emission associated with the LLR component in 
contours overlaid on a velocity integrated map of NH3(3,3) in grey scale. There is a very 
high degree of correspondence between the two maps, confirming that LLR gas traces clouds 
observed in lower rotation inversion transitions. Spectra F-L in Figure 10 show characteristic 
line profiles for LLR gas. The right-hand panel of Figure 12 overlays a velocity integrated 
map of NH3(6,6) emission associated with HLR gas on a velocity integrated map of HCN(1- 
0) from Wright et al. (2001). HLR gas is predominantly located within 2 pc of the Galactic 
center and interior to the CND. A tongue of HLR gas also extends to the east of the nucleus 
(60", —30") with a velocity of 25 km s~^. This HLR gas does not appear to be associated 
with the streamer at 60 km s~^ and offset roughly 50" to the northeast of Sgr A* that was 
detected in HCN(l-O) by Ho (1995). 

Figure 9 shows the location of a position-velocity cut that begins in the southern 
streamer and passes northwards through the central 2 pc of the Galaxy. The corresponding 
position velocity diagram of Si,{6, 6) — 5*1,(3, 3) is shown in grey scale in Figure 13. Although 
S'jy(6, 6) — Su{3, 3) ranges between -0.65 and 0.13, the grey scale stretch ranges from -0.25 to 
0.25 to highlight the different components in the data. Black corresponds to LLR gas, and 
white corresponds to HLR gas. The 0" position corresponds to A6 = 0". 

Panels h — d oi Figure 13 compare the difference map to emission from three molecular 
tracers. In Figure 13b, we overlay a position- velocity cut from our NH3(6,6) data cube. 
Contours are in steps of 3(T66,c/i5 where o"66,c/i = 0.014 Jy Beam~^ (Herrnstein & Ho 2002). 
Emission from NH3(6,6) is present in both the HLR and LLR gas. Figure 13c plots a 
position-velocity cut of NH3(3,3) with contours in steps of 3, 6, 12, 24, and 48cr33,c/i, where 
(733,cft = 0.011 Jy Beam-i (McGary et al. 2001). Unlike NH3(6,6), NH3(3,3) emission is only 
present in the LLR gas. The lack of NH3(3,3) emission is consistent with all of the HLR gas 
resulting from absorption of NH3(3,3) by cool intervening material along the LOS. 

In Figure 13d, we plot a position- velocity diagram of HCN(l-O) from Wright et al. 
(2001). This position-velocity diagram resembles neither the NH3(6,6) nor (3,3) plot. From 
— 150" to —30", only the line wings of the southern streamer are detected. Self- absorption of 
HCN(l-O) is detected throughout the central 10 pc of the Galaxy (see spectra A, B, C, D, 
and F in Wright et al. (2001)) and strongly affects emission from the massive clouds (McGary 
et al. 2001). The strong self-absorption results from the low equivalent temperature of the 
HCN(l-O) transition (5 K). The higher equivalent temperatures of NH3 (equal to 125 K 
for NH3(3,3)) make the rotation inversion transitions much less susceptible to absorption. 
Although not precisely coincident, the HCN(l-O) emission appears to be most closely related 
to emission from NH3(3,3) and the LLR gas. 
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6.1. The Circumnuclecir Disk 

Faint NH3(3,3) emission associated with the CND is detected at —30 km s^^ at —20" in 
Figure 13c. The CND is the dominant feature in the HCN(l-O) data, with emission observed 
from -70 km s"^ at -30" to +100 km s"^ at 80" (see Figure 13d). However, the HCN(l-O) 
emission is not coincident with the NH3(3,3) emission. In particular, HCN(l-O) emission 
at —20" shows two peaks at —20 and —40 km that bracket the NH3(3,3) emission at 
—30 km s~^, indicating that self-absorption of HCN(l-O) is a significant effect even in the 
CND. 

It is well-established that there is a lack of HCN(l-O) emission in the southeastern 
part of the CND (Giisten et al. 1987; Wright et al. 2001). In Figure 13d, we find that 
emission from HCN(l-O) and HLR gas do not overlap at any position. Therefore, it might 

be suspected that the same material that absorbs NH3(3,3) to produce un-physical (6,6)-to- 
(3,3) hne ratios also absorbs HCN(l-O) emission from the eastern side of the CND. However, 
a detailed comparison of the HCN and NH3 data in the "missing" southeastern part of the 
CND indicate that this scenario is unlikely. Both HCN(l-O) and HLR gas are absent at a 
position of (40", -10") relative to Sgr A* (see Figure 12). At (25", -20") HLR gas overlaps 
part of the missing ring, but it has velocities from 30 — 100 km s~^ and will not absorb 
emission from the CND, which is expected to have velocities between ±20 km s~^(see Figure 
12 of Wright et al. (2001)). Finally, lower hne ratios near 20 km s~^ indicate that the HLR 
gas cloud has a small filling factor compared to the southern streamer. It is thus unlikely 
that the cloud could absorb a significant portion of the CND. The "missing" southeastern 
part of the CND most likely reflects an intrinsically clumpy morphology. 



6.2. The Southern Streamer 

The majority of LLR gas in Figure 13 is associated with the southern streamer. This 
northeastern extension of the 20 km s~^ CMC towards the CND has been suggested to be 
physically connected to the nucleus (Ho et al. 1991; Coil & Ho 1999; McGary et al. 2001). 
However, this model has been hindered by a small velocity gradient and only minimal heating. 
For a central mass of 10^ Mq, a cloud at a distance of 5 pc should have a velocity gradient of 
~ 9 km s^^ pc^^ (Coil & Ho 1999). However, a velocity gradient of only 2-3 km s^^ pc^^ is 
measured along the 5 pc length of this cloud (Coil & Ho 1999). The small observed velocity 
gradient imphes that the bulk of the acceleration must occur perpendicular to the LOS. This 
orientation seems unhkely, however, because the 20 km s~^ CMC is known to be located in 
front of the Galactic center along the LOS (see §7). 
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LLR gas associated with the southern streamer has a typical velocity of 20-30 km s~^. 
Strong NH3(3,3) emission is associated with this feature (see Figure 13c), and even the outer 
satellite hyperfine lines, which are offset by ±29 km s~^ from the main line, are detected. 
Emission from the southern streamer appears to continue to the north of the CND, with no 
evidence for a significant change in velocity. LLR gas can also be seen in Spectrum A of 
Figure 10. In this spectrum, the LLR gas shows slightly higher velocities of ~ 50 km s~^. 
This molecular gas at 50 km is most likely responsible for the stimulated H 92a emission 
observed towards Sgr A* by Roberts & Goss (1993). Assuming that the LLR gas in Spectrum 
A is associated with the southern streamer, we measure an upper limit for the velocity 
gradient along the streamer of 12 km s~^ arcmin"^ (5 km s~^ pc~^), which is still a factor 
of two smaller than the expected gradient for infalling gas. Increased line widths near the 
CND were also reported by (Coil & Ho 1999). In NH3(3,3) we see no significant evidence 
for line width broadening associated with the LLR gas. Based on the continuation of the 
streamer to the north of Sgr A* as well as the lack of a large velocity gradient or increased 
line widths, we conclude that the southern streamer shows no conclusive evidence for direct 
interaction with material in the central 2 pc of the Galaxy. Instead, its apparent approach 
towards the Galactic center is most likely the result of projection along the LOS in front of 
Sgr A*. 

Although we favor the above model for the southern streamer, the 15" resolution of our 
data prevents us from unequivocally ruling out an interaction of some type with the nucleus. 
While the kinematics of the southern streamer do not indicate a close association with the 
nucleus, it should be noted that a significant change in the NH3(3,3) emission occurs at 
roughly the location of the southern lobe of the CND. NH3(3,3) emission north of the CND 
is significantly weaker than the emission to the south. The coincidence of this change with 
the nucleus may indicate that there is some interaction between this cloud and the nucleus. 

It is also possible that the faint LLR emission to the north of Sgr A* is not associated 
with the southern streamer at all. LLR emission from — 50" in Figure 13c appears to fill 
in the velocities between the two lobes of high velocity emission associated with the CND. 
Continuous emission from —30" to +80" would imply that the CND extends very close to 
the nucleus, which is contrary to the current model of the CND as a ring-like structure. 
However, observations of HCN(l-O) emission from the CND have shown that it is warped 
and clumpy, and it differs significantly from a simple Keplerian ring (Giisten et al. 1987). 
It has even been suggested that the structure is the result of three independent clouds in 
orbit about the nucleus (Wright et al. 2001). A more precise understanding of the CND is 
necessary to rule out an association with the 30 km s~^ LLR gas to the north of Sgr A*. 

Finally, the faint NH3(3,3) emission to the north of Sgr A* may partially result from 
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"blooming" of flux during the MEM deconvolution. (Blooming occurs when extra flux in the 
form of a faint, extended component is added to the image by MEM in order to obtain the 
"best" solution.) Because emission to the north of the CND has slightly more red-shifted 
velocities than the bright southern emission, we believe that it is unlikely that this emission 
is the result of blooming. These questions should be resolved with the future addition of 
single-dish data, which will eliminate deconvolution artifacts and also allow us to use the 
full 3" resolution of the data. 



6.3. The High-Line-Ratio Cloud 

High line ratio gas has a very different kinematic structure from the southern streamer 
and CND. Figure 13 shows that HLR gas has a large velocity gradient from ^ 50 km s^^ 
at -20" to -20 km s"^ at ^ 0" and then back to km s"^ at -^20". Due to the large 
range of velocities over which HLR gas is observed, a cool gas cloud along the LOS cannot 
produce the majority of the absorption of NH3(3,3). The absorbing material responsible for 
the un-physical (6,6)-to-(3,3) line ratios most likely lies in a shielded outer layer of the same 
cloud. The rotation of the HLR cloud is in the opposite sense as the CND. The gradient 
is in the same direction as the high velocity cloud discussed in McCary et al. (2001), but 
the positions and velocities of the emission from the two features do not precisely agree. 
Improved spatial resolution is necessary to determine the exact morphology of this cloud. 

The large velocity gradient in the HLR gas explains the observed line widths of of 80- 
90 km s~^ in Figure 2 of Herrnstein & Ho (2002). The true line width of the gas can be 
measured at the turning point in the NH3(6,6) position- velocity diagram (position of 0" in 
Figure 13b). We estimate that the intrinsic line width of the HLR gas is ~ 20 km s~^. This 
line width is very similar to intrinsic line widths determined for the other molecular features 
in the central 10 pc (see Figure 2c). 

Line ratios at the position where HLR gas crosses the 20 km s~^ CMC favor a model in 
which the HLR cloud is significantly smaller than the synthesized beam of these data. At a 
position of —15" in Figure 13a, there is a slight decrease in the values of 5',y(6, 6) — 5,^(3, 3). 
Lower (6,6)-to-(3,3) line ratios are consistent with beam dilution of line ratios at these 
velocities by cool emission from the southern streamer. 

The "hook" profile of the position-velocity diagram for the HLR gas is characteristic of 
a cloud that is undergoing both rotation and either infall or expansion (Keto et al. 1988). 
The position-velocity profile of the HLR gas is best fit by a hnear velocity gradient of - 
1.6 km s~^ arcsec"^ (—40 km s~^ pc~^) and expansion or contraction at 80 km s"^. Assuming 
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the gradient is due to circular rotation and that the measured length of the cloud (~ 1 pc) is 
twice the radius, the observed gradient corresponds to rotation at 20 km s~^. This rotation 
velocity implies a minimum enclosed mass of 4 x 10^ Mq. If the cloud is assumed to have a 
circular orbit around a mass of 10'' Mq, then the observed velocity gradient would correspond 
to an inclination of ~ 4°. Such a low inclination is unlikely, although it cannot be ruled out 
with the current resolution, and the HLR cloud is probably not in a circular orbit. 

Because emission is blue-shifted from the expected rotation velocities, we must either be 
observing emission from the far side of a contracting cloud or the near side of an expanding 
shell. We currently favor a model in which we arc observing the near edge of an expanding 
shell. Assuming the source of the expansion is therefore behind the observed shell along 
the line-of-sight, then it naturally follows that the cool layer which absorbs the NH3(3,3) 
would be located between the hot material on the inner edge of the shell and the observer. 
Although it is possible that the other half of the shell is truly missing, the lack of symmetry 
may also provide clues to the morphology of the shell and the location of the cool absorbing 
material. Improved spatial resolution will be useful in resolving this question. 



7. Towards a 3-D Model of the Central 10 pc 

In the preceding sections, we have investigated the physical conditions and interactions 
of dense molecular clouds within 5 pc of Sgr A*. We briefly summarize our main results 
below: 

• The 50 km s^^ GMC shows no evidence for a strong interaction with Sgr A East. 
Based on the observed effects of Sgr A East on the western streamer, we find that the 
energy of Sgr A East is consistent with a single supernova occurring roughly 10"^ years 
ago (§4). 

• The temperature structure of dense clouds in the central 10 pc can be characterized 
by a two-temperature structure on scales of 0.5 pc. The data are consistent with the 
structure on scales of 3 pc proposed by (Hiittemeister et al. 1993a), in which roughly 
one quarter of the gas is in a hot component at 200 K. 

• NH3 (6,6)-to-(3,3) line ratios in the central 2 pc exceed theoretical limits for gas in 
LTE (§5.2). We conclude that these line ratios result from absorption of NH3(3,3) by 
a cool, shielded layer of the same cloud. This result implies that n^/j ~ 10^ cm~^. 

• Three distinct molecular gas components lie within a projected distance of 2 pc from 
Sgr A* (§6). Both the CND and the "high-line-ratio gas" show large velocity gradi- 
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ents, although in opposite directions. Emission associated with the southern streamer, 
however, shows no significant velocity gradient, and extends to the north of the nucleus. 

It is important to relate the molecular features described in this paper to other structures 
in the nuclear region. Discussions of the central parsecs of the Galaxy seem to inevitably 
evolve into discussions of the location of these structures along the LOS. In this concluding 
section, we present a model for the LOS arrangement of the features in the central 10 pc. The 
model is based on the molecular line data presented in this paper as well as the most current 
results in the published hterature. The model is not to scale (or completely inclusive), but 
it is presented to promote discussion of the interaction between different structures observed 
at the Galactic center. 

A schematic drawing of the main structures in the central 10 pc of the Galaxy is pre- 
sented in Figure 14. The elongated Sgr A East shell fills the center of the image. Towards 
the western side of the shell, a black dot shows the location of Sgr A*, which is surrounded 
by the arcs of Sgr A West (light grey) and the backwards "C" shape of the CND (black). 
The cloud of HLR gas discussed in §6 is also shown in close proximity to Sgr A*. Because its 
exact morphology is not known, it is shown as a grey circle, slightly offset to the southeast 
of Sgr A*. The major molecular features described in this paper are labeled. Hatched lines 
extending northwards from the southern streamer indicate the faint extension of this cloud 
to the north of Sgr A* (see §6.2). A second supernova, G359.92-0.09, which is impacting 
Sgr A East in the southeast is also shown (Coil & Ho 2000). Positions of 1720 MHz OH 
masers (Yusef-Zadeh et al. 1999) are marked with crosses, and the Hll regions observable at 
6 cm are shown as small circles in the 50 km s~^ GMG. 

In Figures 15 and 16, we present schematic drawings of the proposed location of these 
structures along the LOS. Figure 16 also shows the motions of the clouds inferred from their 
Doppler-shifted velocities. The LOS structure of the arcs of Sgr A West is based on H 92q; 
emission from Roberts & Goss (1993). Because the location of the eastern bar relative to 
Sgr A* is not well known, it is not included in the LOS model. The northern arc is best 
described as a tidally disrupted cloud that is passing within 0.13 pc of Sgr A*. The orbit 
produces high negative velocities of -250 km s~^ at the southern-most end of the arc (Roberts 
& Goss 1993; Roberts et al. 1996). The western arc, however, is coincident with the inner 
edge of the circumnuclear disk, and most likely is the ionized inner layer of this structure 
(Roberts & Goss 1993). Absorption of starlight and Br7 emission from the western arc of 
Sgr A West favors an orientation of the CND with the western side in front of the Galactic 
center along the LOS (Roberts & Goss 1993). The location of the HLR molecular gas is 
schematically drawn interior to the CND, but exact determination of the morphology and 
location of this cloud cannot be determined with the current resolution of our data. 
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Absorption of 90 cm continuum emission from Sgr A East by Sgr A West places the 
nucleus in front or just inside Sgr A East (Pedlar et al. 1989). The connection of the northern 
ridge, which lies along the northern edge of Sgr A East, to the CND indicates that Sgr A 
East must be within a few parsecs of the nucleus (McGary et al. 2001). More recently, 
Macda ct al. (2002) place Sgr A*, Sgr A West, and the CND just inside the leading edge of 
Sgr A East. (For schematic drawings of proposed relative locations of Sgr A West, Sgr A*, 
and Sgr A East, see e.g. Pedlar et al. (1989) and Maeda et al. (2002).) 

As discussed in §4, the velocity gradient of the western streamer indicates that it is 
highly inclined to the LOS and expanding outwards with Sgr A East. The northern ridge is 
placed along the northern edge of Sgr A East. The cloud is located at the same distance as 
the center of Sgr A East, with an orientation roughly in the plane of the sky. This placement 
allows for the bulk of the motion associated with the expansion of Sgr A East to occur 
perpendicular to the line of sight, consistent with its velocity of -10 km s~^. 

The 50 km s~^ cloud is located predominately to the east of Sgr A East. As discussed 
in §6, stimulated H 92q; emission at 50 km s~^ observed towards Sgr A* (Roberts & Goss 
1993) is most likely associated with the southern streamer and not the GMC. The 50 km s~^ 
cloud is not observed to strongly absorb X-ray emission from the central 10 pc (Park et al. 
2004), which also implies that little of the material is in front of the Galactic center along 
the LOS. This GMC is connected to the 20 km s~^ cloud along the molecular ridge (Coil & 
Ho 2000). As mentioned in §4, the Hll regions found along the interior of the western edge 
of this GMC are much too old to have been caused by the current impact of Sgr A East. 

Observations of formaldehyde absorption towards the Galactic center, indicate that the 
20 km s^^ GMC lies in front of the nucleus along the LOS (Giisten & Downes 1980). This 
placement has been recently confirmed by absorption of 2-10 keV X-rays in the region of 
this GMC (Park et al. 2004). The distance between the 20 km s^^ GMC and the nucleus is 
constrained by observations of SNR 359.92-0.09. This SNR was first reported by Coil & Ho 
(2000), who detected it in 20 cm radio continuum images from Yusef-Zadeh & Morris (1987). 
More recently, an X-ray filament associated with the southwest edge of this SNR has been 
detected with both XMM-Newton and Chandra (Sakano et al. 2003; Park et al. 2004). The 
impact of SNR 359.92-0.09 is responsible for the concave morphology of Sgr A East in the 
southeast. Multiple 1720 MHz OH masers with velocities of ~ 50 km s~^ have been detected 
along this interface region (Yusef-Zadeh et al. 1999). These masers have been associated 
with disrupted gas associated with the 20 km s~^ GMC (Coil & Ho 2000). Because SNR 
359.92-0.09 is interacting with both Sgr A East and the 20 km s"^ GMC, the GMC must 
be no more than 5 — 10 pc from Sgr A*. Finally, based on the results of §6, the southern 
streamer is depicted as extending northwards from the 20 km s"^ cloud, but not interacting 
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with the nucleus. 

Although uncertainties in this model undoubtedly remain, the fact that such a detailed 
three dimensional model of the Galactic center can be proposed based on observational data 
reflects the great advances in recent years in our understanding of this region. In particular, 
the data presented in this paper have shed new light on the role of dense molecular gas in 
the region around the supermassive black hole. The large velocity coverage of our VLA data 
combined with the relatively high excitation requirements of the rotation inversion transitions 
have resulted in a more complete picture of the location of dense material within 5 pc of 
Sgr A*. The additional knowledge of temperatures, line widths, and cloud masses, all of 
which are directly obtained from the NH3 data, have enabled us to more precisely determine 
the interactions between these dense molecular clouds and other features near the nucleus. 
At present, a more precise understanding of the role of the dense molecular structures in the 
Galactic center is primarily hindered by the relatively low resolution (16" x 14") of our data. 
Future incorporation of measurements made with large single-dish telescopes should allow 
us to resolve many of the outstanding questions regarding dense molecular gas at the center 
of the Galaxy. 

The authors would like to thank R. Zylka for the 1.2 mm continuum image, M. Wright 
for the HCN(l-O) data, and the referee for a very careful and thorough report. R. M. H. 
would also like to thank J. Greene and P. Sollins for many productive discussions. The 6 cm 
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Tablc 1. Physical Parameters for Galactic Center Clouds' 



Parameter 


Unit 


Southern 


Northern 


SEl 


SE2 


Northern 


50 km s-i 


Western 






Streamer 


Fork 






Ridge 


GMC 


Streamer 


Extent in R.A. 


// 


-30, 60 


-70, -30 


60, 90 


90, 140 


30, 100 


80, 150 


-100, -40 


Extent in Dec. 


// 


-30, -150 


-150, -90 


-30, -140 


-10, -100 


80, 130 


0, 125 


-100, 100 


(A Vint) 


km 


15.2 


13.9 


14.6 


16.7 


15.8 


16.9 


20.5 


O-(Avint) 


km 


8.1 


5.4 


8.1 


10.0 


4.0 


6.4 


4.5 


(T-m(l,l)> 




2.4 


3.0 


2.2 


1.9 


0.5« 


2.5 


0.4« 


<^(Tm(l,l)) 




0.7 


0.2 


0.7 


0.3 


0.4" 


0.5 


0.3" 


(T«21> 


K 


26.4 


25.4 


28.9 


29.3 


28.2 


21.4 


30.9 


(r{TR2l) 


K 


5.1 


2.8 


5.2 


4.7 


4.5 


2.9 


3.9 


{Tk)'° 


K 


33 


32 


39 


40 


37 


25 


46 


{Tex) 


K 


6.2 


4.3 


4.9 


4.7 


4.9'' 


4.4 


4.1" 




K 


1.1 


0.2 


0.5 


0.6 


lA" 


0.5 


0.7" 


(^(1.1) > 


10^^ cm-2 


3.5 


3.6 


1.5 


1.2 


0.9 


2.8 


1.4" 


'^(^^(1,1)) 


IQlS cm-2 


1.9 


1.0 


1.1 


0.9 


0.2 


1.1 


0.6" 




IQlS cm-2 


8.7 


5.0 


3.0 


2.9 


2.1 


6.1 


3.7" 




10^'' cm"2 


4.3 


4.2 


2.1 


2.1 


0.5 


2.6 


1.4" 




lo'^n" 


4.0 


n.3 


1.5 


2.2 


0.4 


3.3 


0.4 


Mass'' 


..1/ . 


<s X 10' 


i X lO'' 


ixlO' 


2 X iO ' 


2 X iO-'' 


X 11)-^ 


i X 10'^ 



^The mean and standard deviation of each parameter is calculated using all pixels where a value was determined from 
Monte-Carlo simulations with > 3cr significance. 

''(Tff ) is calculated from (Tfl2i) using Danby et al. (1988), which assumes n(H2) = 10^ cm~*. 

'^Cloud area is calculated from the number of pixels with Nj^jhs determinations (see Figure 5a). 

"^Masses are calculated assuming X(NH3)= lO"'^. Because the value of X(NH3) is not well-determined, masses should be 
taken as order of magnitude estimates. 

"No pixels with > 3a significance for this feature. All pixels with a determination from the Monte-Carlo simulation are 
included in the calculation. 
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Table 2. Summary of NH3 Observations at the Galactic Center 



Label 


Telescope 


Beamt (") 


Measured Temperatures 


Reference 


G81 


Effelsberg 


40 


Tr21, Tr4i, Tr42 


Glisten et al. (1981) 


089 


NMA 


24 X 13 


Tr21 


Okumura et al. (1989) 


ASS 


Haystack 


84 


Tr21, Trqs 


Armstrong & Barrett (1985) 


H93 


NRAO 40-m 


84 


Tr21, Tr54 


Hiittemeister et al. (1993a) 


C99/00 


VLA 


14 X 9* 


Tr21 


Coil & Ho (1999, 2000) 


HOI 


VLA 


15 X 13 


T 1-121 


this ])a.])er 



^For interferometric observations, this synthesized beam is listed. 

^Coil & Ho (1999, 2000) report average values for each condensation, which is typically 
30" in size. 
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Table 3. NH3 Rotational Temperature Measurements for Galactic Center Molecular 

Clouds 



Aa 


AS 


Tr21 


Tr41 


Tr42 


Tr54 


Tr63 


Ref. 


(") 


n 


(K) 


(K) 


(K) 


(K) 


(K) 




ioo 






> 4/ 










76 


95 




> 44 


> 41 






G81 


143 


53 




> 52 


> 54 






G81 


152 


47 




> 49 


> 56 






G81 


177 


32 




> 39 


> 43 






G81 


122 


19 




> 49 


> 61 






G81 


62 


-155 


66'!" 


46 


46 






G81 


30 


-207 






49 






G81 


-3 


-259 


o5_io 


451}^ 


49lg 






G81 


-67 


-287 


96 


47 


43 






G81 


-33 


-308 


lOOt^^ 




37l^ 






G81 


-66 


-360 


71 


46 


43 






G81 


-13 


-392 


63 


44 


46 






G81 


-97 


-412 


91 


40 


37 






G81 


39 


-424 


51 


41 


46 






G81 


219 


100 










80 ± 16 


A85 


94 


70 










80 ± 16 


A85 


172 


23 


34 ± 7 








75 ± 15 


A85 


125 


-54 


37 ± 7 








75 ± 15 


A85 


77 


-130 


36 ± 7 








100 ± 20 


A85 


30 


-207 


48 ± 10 








80 ± 16 


A85 


-18 


-283 


31 ± 6 








70 ± 14 


A85 


59 


-331 










100 ± 20 


A85 


-66 


-360 


62 ± 8 








70 ± 14 


A85 


-113 


-436 










90 ± 18 


A85 


-8 


-203 


25 ± 7 










089 


-38 


-203 


24 ± 5 










089 


-18 


-233 


31 ± 7 










089 


-49 


-253 


30 ± 13 










089 


-77 


-253 


22 ± 25 










089 


-18 


-263 


24 ± 3 










089 


-28 


-263 


17 ± 2 










089 
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Table 3 — Continued 



Aa 


AS 


Tr21 




Tr42 


7r54 


^^63 


Ref. 


n 


n 


(K) 


(K) 


(K) 


(K) 


(K) 




-70 


-263 


27 ± 34 










089 


-58 


-323 


18 ± 3 










089 


-107 


-323 


3 ± 6 










089 


-91 


-347 


18 ± 3 










089 


-49 


-353 


13 ± 3 










089 


-106 


-335 


35 ± 2 






64 ±4 




H93 


26 


-23 


32+f 










C99 


26 


-83 


22+f 










C99 


-13 


-163 


22+f 










C99 


144 


4 


20+f 










COO 


105 


-126 


30+f 










COO 



'''Rotational temperature measurements by Giisten et al. (1981) 
do not include an estimate of uncertainties. For these measure- 
ments, we assume an uncertainty of ±10 K. 
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Fig. 1.— Velocity integrated maps of a) NH3(1,1), b) NH3(2,2), c) NH3(3,3), and d) 
NH3(6,6). Contours are in steps of 1.32 Jy beain"^ km s~^ in all four panels. The rms noise 
level, ajK, for each (J, K) rotation inversion transition is an — 0.28 Jy beam"-*^ km s~^, (T22 = 
0.30 Jy beam~^ km s~^, (733 = 0.33 Jy beam~^ km s~^, and a^e — 0.33 Jy beam"^ km s"^ 
(McGary et al. 2001; Herrnstein & Ho 2002). The position of Sgr A* is labeled by a star in 
each panel. The main molecular features in the central 10 pc are labeled in panel c. 
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Fig. 2. — NOTE: Image too large for astro-ph. Please download a full version of this paper at 
http://www. astro. Columbia. edu/^herrnstein/NH3/paper/ .a) Mean value of NH3(1,1) main 
line opacity (Tin(l, 1)). Contours are in integer steps, b) Significance of opacity determina- 
tions. Contours are 1, 2, 3, 6, 9, and 12cr. Hatches denote pixels for which only a lower 
limit for rin(l, 1) could be determined, c) Mean value of intrinsic line width (Avint)- Typical 
intrinsic line widths are 15 km s~^. d) Significance of intrinsic line width determinations. 
Contours are in steps of 3a. Hatches denote pixels for which only an upper limit for Avint 
could be determined. 



Fig. 3. — NOTE: Image too large for astro-ph. Please download a full version of this paper 
at http:/ /www. astro. Columbia. edu/'^herrnstein/NHS /paper/ . a) Pseudo-color map of the 
derived (2,2)-to-(l,l) rotational temperature. Contours are in steps of 10 K. b) Significance 
of Tr2i determinations. Contours are in steps of 3(T. c) Derived values for the excitation 
temperature. Contours are in steps of 4 K. Compared to the rotational temperature, the 
excitation temperature is smoothly distributed (especially in regions with high S/N) with 
values between 4 and 8 K. d) Significance of the T^^ values. Contours are 3, 15, 30, and 45cr. 
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103 ^QA 106 10^ 

□ (H^) [cm"'] 

Fig. 4. — Excitation temperature as a function of H2 density for kinetic temperatures of 10, 
30, 100, and 300 K. This plot is valid for NH3(1,1), (2,2), (3,3), and (6,6) because the critical 
densities for these transitions are all equal to 2 x 10^ cm~^. Gas is not thermalized until 
n{H2) ^ 10^ cm-^. 



Fig. 5. — NOTE: Image too large for astro-ph. Please download a full version of this paper 
at http: //www. astro. Columbia. edu/^herrnstein/NH3/paper/ . Total column density of NH3. 
Contours are in steps of 4 x 10^^ cm~^. Significance of A^nhs determinations. Contours are 
in steps of l.Scr. 
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Fig. 6. — Calculated kinetic temperature for Tr2i, 7k42; ^_R54; and Tr63 ^or six different 
combinations of hot gas at 200 K and cold gas at 15 K. Weighted mean kinetic tempera- 
tures calculated from the Tr21 and T^es measurements presented in this paper are overlaid 
as squares, b) Weighted mean kinetic temperatures calculated from published Tr2i, 2^_r42, 
and Tr63 measurements in Figure 8 are overlaid as filled circles. The measurements of 7^21 
and Tij54 from Hiittemeister et al. (1993a) are overlaid as triangles. 



Fig. 7. — NOTE: Image too large for astro-ph. Please download a full version of this paper 
at http: //www. astro. Columbia. edu/^herrnstein/NH3 /paper/ . a) Position of the single dish 
measurements listed in Table 3 overlaid on a grey scale image of 1.2 mm emission (Zylka 
et al. 1998). Contours denote 10, 20, 40, and 80% of the peak flux. Labels correspond to the 
list of observations in Table 2. b) Position of interferometric observations in the same region. 
Primary beams of the interferometer are plotted. For Okumura et al. (1989), the synthesized 
beam of the NMA is also overlaid at the positions where temperature measurements were 
made. 
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Fig. 8. — Measured kinetic temperature versus distance from Sgr A* calculated from the 
rotational temperatures listed in Table 3 using Danby et al. (1988). The four panels show- 
kinetic temperatures derived from a) Tji2i, b) Tr42, c) Tf^r^i, and d) T/^es- Uncertainties in 
the radial distance from Sgr A*, which are equal to the FWHM of the beam, have been 
omitted for clarity. The dark grey line in panel a shows the best fit radial profile for kinetic 
temperatures derived from our data. The surrounding light grey box shows the approximate 
range of kinetic temperatures from pixels with a significance of > 3a. 
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Fig. 9. — Position of twelve spectra and one position-velocity cut overlaid on a velocity 
integrated map of NH3(6,6). Sgr A* is located at (Aa = 0",A5 = 0") (spectrum A). 
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Fig. 10. — Spectra from NH3(3,3) (grey) and NH3(6,6) (black) at the positions labeled 
in Figure 9. For both transitions, we use image cubes deconvolved using a 14" Gaussian 
taper. Spectra are Manning smoothed and the noise levels are cr^^ ch ~ "^-^ Beam~^ and 
(^m,ch = 9-5 mJy Beam-\ 
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Fig. 11. — Schematic model of the orientation of a hot and cold cloud that can produce 
"un-physical" NH3 (6,6)-to-(3,3) hne ratios. 
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Fig. 12. — Left: Velocity integrated map of NH3(6,6) emission at all points with NH3 (6,6)- 
to-(3,3) line ratios less than one overlaid on a grey scale map of velocity integrated NH3(3,3) 
emission. Contours are in steps of 3(766 where — 0.33 Jy Beam~^ km s~^. Low line ratio 
gas has a similar distribution to velocity integrated emission from lower-energy rotation 
inversion transitions. Right: Velocity integrated map of NH3(6,6) emission at all points 
with NH3 (6,6)-to-(3,3) line ratios greater than one overlaid on a grey scale map of velocity 
integrated HCN(l-O) emission (Wright et al. 2001). Contour levels are the same as in the 
left-hand panel. High-line-ratio (HLR) gas is confined to the inner 2 pc, interior to the CND, 
and a tongue of emission extending to the east at A6 = —30". 
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Fig. 13. — Position velocity diagrams showing the distribution of molecular gas in the central 
2 pc. The position of the cut is labeled in Figure 9. A position of 0" corresponds to AS = 0", 
and positive positions are to the northwest of Sgr A*. In all four panels, Siy{6, 6) — 5*1,(3, 3) 
is plotted in grey scale. High- line-ratio gas (5',^ (6, 6) > Si,{3, 3)) appears as white, while low- 
line-ratio gas {Si,{6, 6) < S^{3, 3)) is dark grey. Although Si,{6, 6) — Si,{3, 3) ranges between 
-0.65 and 0.13, the grey scale stretch ranges from -0.25 to 0.25 to highlight the different 
components in the data, b) NH3(6,6) overlaid in contour steps of 3(T66,c/i- NH3(6,6) emission 
is observed from both the HLR and LLR gas. c) NH3(3,3) overlaid in contour steps of 3, 
6, 12, 24, and 48 cr33,ch- NH3(3,3) emission comes only from LLR gas. d) HCN(l-O) from 
Wright et al. (2001) overlaid in contour steps of 0.25 Jy Beam~^. The lack of HCN(l-O) 
emission at 30 km s~^ is the result of self-absorption. 




Fig. 14. — Schematic drawing of the Galactic center as seen in the plane of the sky. 
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Fig. 15. — Schematic drawing of the Galactic center with left denoting the direction to the 
observer and the vertical direction denoting north and south. 
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Fig. 16. — Schematic drawing of the Galactic center with left denoting the direction to the 
observer and the vertical direction denoting east and west. Arrows show the inferred motions 
of the main features based on the Doppler-shifted velocities. 



